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Abstract 



We use the QCD sum rule approach to calculate the splitting between vector 
and pseudoscalar mesons containing one light and one heavy quark, and the 
kinetic energy of the heavy quark. Our result for the splitting induced by the 
chromomagnetic interaction agrees to the experimental data on charm and 
beauty mesons. For the matrix element of the kinetic energy operator, we 
obtain the value K = -(0.60 ± 0.10) GeV 2 . 
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I. INTRODUCTION 



In recent years there has been a continuous interest in the study of mesons built of one 
heavy and one light quark. This interest is fuelled by a constant flow of new experimental 
data and the results of lattice calculations. The main theoretical achievement in the past few 
years was the development of the heavy quark effective theory (HQET, see |J for reviews), 
which provides one with a systematic tool for the study of the heavy quark limit rriq — > oo 
and for the classification of corrections which are suppressed by powers of the heavy quark 
mass. Physical applications of this approach include various decays of charm and beauty 
mesons, for which the preasymptotical corrections ~ l/mg are likely to be significant. 
Considerable effort has been made in the past two years to estimate them, see e.g. 
but this task is far from being completed yet. The Lagrangian density, written in terms of 
the effective heavy quark fields, contains to 0(l/rriQ) accuracy two additional contributions 
|||7J apart from the leading one. They are related in an obvious way to the non-relativistic 
kinetic energy operator, and to the Pauli term, describing the chromomagnetic interaction. 
Matrix elements of these two operators over meson states are fundamental observables in 
the heavy quark effective theory, and determine the next-to-leading order corrections to 
the meson masses suppressed by powers of the quark mass. The matrix element of the 
chromomagnetic interaction operator is a leading contribution to the mass splitting between 
vector and pseudoscalar mesons, and can be directly related to the data on the meson 
spectrum. The matrix element of the kinetic energy operator contains information about 
the smearing in heavy quark momentum, which is important in many applications and 
contributes, e.g., to subleading form factors of semileptonic decays ||, and recently was 
shown to give a significant contribution to the electron spectrum in inclusive B-decays |J . 

Up to now, there exist no quantitative estimates for the kinetic energy of the heavy 
quark, except for a single attempt in M. This task proves to be difficult for lattice QCD, 



because of power divergences which need to be subtracted, see e.g. [10]. In this paper we 



calculate the kinetic energy and the chromomagnetic mass splitting using the QCD sum rule 



approach [pTlJI , including radiative corrections. 

The presentation is organized as follows. Sec. II is introductory, and contains a short 
discussion of the heavy quark effective theory, giving necessary definitions and notations. 
Some results of more technical nature are given in App. |A]. The QCD sum rules for the 
relevant three-point functions in the effective theory are derived in Sec. III. Sec. IV contains 
our main result, namely the sum rule for the matrix element of the kinetic energy operator, 
including two-loop radiative corrections and renormalization group improvement to two- 
loop accuracy. Details of the calculation are given in Apps. and 0. Finally, in Sec. V we 
give the summary and conclusions. 



II. THE HEAVY QUARK EXPANSION 

The dynamics of hadrons containing both light and heavy quarks can be described by 
an effective field theory, in which the heavy degrees of freedom are integrated out, and 
the resulting effective Lagrangian is expanded in inverse powers of the heavy quark mass. 



Following jl2| we introduce the heavy quark effective field h v as 
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P+Q(x) = e- im ^ v - x) h v (x) } (2.1) 

where 

P + = V) (2-2) 

is the projector on the upper components of the heavy quark field Q(x), vtlq is the heavy 
quark mass, and v its four-velocity. The contributions of lower components of the Q-quark 
field are integrated out, producing an effective theory with the Lagrangian density 

£ = h v i{v ■ D)h v + —L/C + — L <S + 0(l/m%), (2.3) 
IrrtQ ZrriQ 

where 

K = KiiD^fK (2.4) 
is the operator of the non-relativistic kinetic energy, and 

S 4(w) 3/&W9F ^' (2 ' 5) 

is the Pauli term. The covariant derivative is defined as _D M = <9 M — igA^. In (|2.4j) we have 
introduced projectors on the directions orthogonal to the velocity v^. 

9 hi/ = 9\i.v v^v u , 

(D^f = DpD» - {vDf. (2.6) 

The operator of the chromomagnetic interaction has a non-trivial anomalous dimension, 
and the coefficient in front of it receives corrections of higher orders in the coupling constant, 
Cp{fi = ttiq) = 1/2 + 0(a s (mg)). In contrast to this, the kinetic energy operator has zero 
anomalous dimension to all orders, and its coefficient is exactly one, cf. [13]. This property 
follows from a residual symmetry of the effective theory under non-relativistic boosts, and 
is related to the fact that the kinetic energy term in the non-relativistic Hamiltonian is 
invariant under these transformations. 

Operators in the full theory (QCD) are expanded in a series of operators in the effective 
theory as 

qcd = J2Ci(m Q /fi,a s ({i))0? QET (fi). (2.7) 

i 

The coefficients in this expansion depend on the scale /i according to the renormalization 
group equations, 



C(m Q //i, a s (fi)) = (7(1, a s {m Q )) exp 



s( m Q) 7(0) , ' 
dg 



(2.8) 



where 7(g) is the anomalous dimension of the corresponding operator in the effective theory, 
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7 = 7o 



9 



167T 2 

and (3(g) is the Gell-Mann-Low function 

,2 



+ 7i 



16tt 2 



(2.9) 



P = ~9 



0o 



9' 



167P 



+ Pi 



16tt 2 



(2.10) 



At the scale fi = uiq, the coefficient functions C, are determined from the condition that the 
matrix elements of effective operators, times the appropriate coefficient functions, equal the 
matrix elements of the corresponding QCD operators at this scale to the required accuracy. 
For the practically important cases of vector and axial currents built of one heavy and one 
light quark, the heavy quark expansion reads @,Q 



i -a 



7X I I 71 



2m Q \2mg tt J 



IX 2 71 



2m Q p \[ 



1 a 



(2.11) 



2m Q 7r J 

where Cf = {N^ — 1)/ (2N C ). Invariant functions and matrix elements of operators in HQET 



are most conveniently defined using the so-called trace formalism |L5] , which makes the spin 
symmetries explicit. Following H, we define 



(0\qTh v \M(v)) = -F(ji)Tr{rM(v)}, 



where A4(v) is the spin wave function 



A A I \ i r> f ~^75 f° r J P = , 



(2.12) 



(2.13) 



The relation in (|2.12 ) is valid for an arbitrary structure of Dirac matrices T. To two-loop 
accuracy we have 



\ Olsiji) j 



7o/(2/? ) r 



! Jo hi Pi\ a s (m Q ) -a a (/x) 
SPo Uo Pol vr 



(2.14) 



with the anomalous dimensions flIl],|i~6|,|T7 



70 = -4 , 7i = -254/9 - 56tt 2 /27 + 20n f /9 
A) = ll-2n / /3, Pi = 102-38n//3. 



(2.15) 
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The effective coupling F(fi) depends on the scale in a similar way. Defining physical lepton 
decay constants, e.g. for B-mesons, by 

(0 1 <nV75&| B(q)} = if B q» 

(0\q llx b\B*(q, e)> = f B *m B ^ e ■ q = 0, e 2 = -1, (2.16) 

and taking into account the matching conditions in ( |2.11[ ), one obtains the relations of the 
physical couplings to the effective ones PJl4}|: 



fsV^ = (l ~ 2 -^^i FQi = m B ) + 0(l/m B ), 

f B *V^= U-~^^F(p = m B )+0{l/m B ). (2.17) 

Furthermore, we define matrix elements of the operators of kinetic energy and chromomag- 
netic interaction over heavy-light mesons as 

(M(v)\K\M(v)) = -KTr{M(v)M(y)}, 

(M(v)\S\M(v)} = -d M ^Tr{M(v)M(v)} . (2.18) 
Here du is defined by 

P + a aP M(v)a al3 = 2d M M(v), (2.19) 

yielding 

Note that the normalization of matrix elements is 

(M(v)\h v h v \M(v)} = 2 ■ (-1) Tr{M(v)M(v)} = 2m M , (2.21) 

so that 

(M(v)\K\M(v)) 



K 



(M(v)\h v h v \M(v)Y 
1 (M(v)\S\M(v)) 



d M (M(v)\h v h v \M(v)} 
To 1/toq accuracy the meson masses are given by 



(2.22) 



m M = m Q + A- 7 ^—[K + d M E] i (2.23) 
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where A determines the difference between the quark and meson masses in the heavy quark 
limit, and is a fundamental observable in the heavy quark effective theory. It has received 
a lot of attention recently, see e.g. [Q,|J. The chromomagnetic interaction determines the 
splitting between pseudoscalar and vector mesons: 

m \ - m 2 p = 4S + 0(l/m Q ). (2.24) 



The relation ( |2.23| ) can be derived in a variety of ways and is widely known One 



possibility is, e.g., to consider the 1/rriQ expansion of suitable two-point correlation functions 
in QCD near the particle-type singularity, corresponding to the lowest-lying meson state, 
cf. App. 0. 



III. THE SUM RULES 

The aim of this paper is to obtain quantitative estimates of the matrix elements K and 
S. We use the method of QCD sum rules, see the book which we apply directly to 
suitable correlation functions in HQET. In particular, we consider the following three-point 
correlation functions at zero-recoil: 

i 2 J dx J ^e^-^^^oirlg^ri^^/c^)^^)^^)}^) = 

= -TT{T 1 P + r 2 }T K (u,u'), (3.1) 

i 2 j dx J d|/e-^)--'^)(0|T{g(x)r 1 ^( ; r) < S(0)^( 2 /)r 2g ( 2 /)}|0) = 

= -d M Ti{T 1 P + T 2 }T^(uj,uj'). (3.2) 

Saturating the three-point functions with hadron states, one can isolate the contribution 
of interest as the one having poles in both the variables to and u' at the value to — u' — A: 

E(/i)F 2 0u) 
Ts( "'" ) = 4(A-^)(A-^) + '" 

T ^-^ 4(A™-^) + --- ^ 

where the scale /i is the normalization point of the currents, and F 2 ([i) is the coupling 
(squared) of the effective current to the lowest-lying meson state, defined in ( |2.12j ). This 
coupling, in turn, can be obtained as the residue of the pole at u> = A in the two-point 
correlation function 

i J <i 4 ie"<" I »(o|T{5(i)r 1 / l „(i)/ l „(o)r 29 (o)} |o) = -i'R{rii 5 +r 2 }nH, (3.4) 

n ^ = W^y> + - (3 ' 5) 
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To derive Eqs. ( fT3| ) and (|3.5|) , the following relation is useful [||: 

J2 Tr^MivftTriMiv)^} = -2m M Tr{r 1 P+r 2 }. (3.6) 

polarizations 

For example, for the lowest-lying meson contribution to the correlation function H(uj), one 
obtains 

which is the result given in ( p.5|) . 

The correlation functions H(uj), T^(u, a/), and Tk(uj, uj') can be calculated in the Euclid- 
ian region, for negative u,u', and receive contributions from perturbation theory and from 
vacuum condensates. The results can be written in form of (double) dispersion relations, 




where the spectral densities are subject to direct calculation in HQET: 

p = p P crt (!) + p ® (qq) + p W {^G 2 ^ + p (5) (qagGq) + ... (3.9) 

The relevant Feynman diagrams, which contribute to the three-point correlation func- 
tions to first order in a s , are shown in Fig. 1 for the chromomagnetic operator, and in 
Figs. 2, 3, and 4 for the kinetic energy (up to dimension 4). The graphs missing there turn 
out to be vanishing. Note that the leading contributions to the correlation function Tg, 
Fig. [I], are of 0(a s ). On the contrary, the leading perturbative contribution to the kinetic 
energy is of order 0(1), see Fig. 2(a). 

In this section, we consider the sum rule for the chromomagnetic mass splitting, taking 
into account the set of diagrams in Fig. 1, and the leading order sum rule for the kinetic 
energy, taking into account the perturbative contribution in Fig. 2(a) and the leading non- 
perturbative correction, which in this case is due to the mixed quark-gluon condensate. The 
contribution to the kinetic energy of the quark condensate is of 0(a s ). The full sum rule 
for the kinetic energy, with account for all the contributions in Figs. 2, 3, and 4 and with 
renormalization group improvement to two-loop accuracy, is considered in detail in the next 
section. 

A straightforward calculation yields: 

Pu(s) = ^s 2 - ^(qq)S(s) + ^(qagGq)5"(s) + 0(a s ), 

3 3 
p K (s, s') = -—s 4 5(s - s') + —(qagGq)S(s)S(s'), 
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p s ( s , s ') = ^§ {s ,2 (s - s73)0(s - s') + s 2 (s' - s/3)6(s' - s)} 

- ^(qq){sS(s') + s'5(s)} + 1 (^G 2 ) <5(s - s') - ±(q*gG q }5(s)5( S >). (3.10) 

Following the usual strategy of the QCD sum rule approach, we apply a Borel transformation 
to the correlation functions in both the external momenta uo and uj'\ 

B w (t) T ^— = lexp(-A/t). (3.11) 
A — oj t 

The transformation B u (t) introduces the Borel parameter t instead of the external momenta. 

We next equate the representations of the correlation functions in terms of hadronic 
states, Eqs. ( |3.3|) and ( |3.5| ), and in terms of spectral densities, Eqs. (|3.8|) and ( |3 . 1 0|) . In the 



latter we constrain the region of integration to the interval of duality < s, s' < Due to 
the symmetry of the correlation functions, it is natural to take the two Borel parameters in 
the three-point functions equal to each other, t% = t<i- The normalization of the Isgur-Wise 
function at zero recoil |19[], further requires that the value of the Borel parameter in the 



three-point function takes twice the value as in the two-point function, t\ — ti — 2t, and 
that the values of the continuum threshold u>o in the sum rules for two- and three-point 
functions coincide. We end up with the set of sum rules[] 



1 ^0 



o 

WO Wo 





wo wo 



^F 2 (/i)S(/i)e- A /* = jdsjds'e-( s+s 'VWp s (s,s'). (3.12) 







The normalization scale \x should be taken to be of order of the typical Borel parameter. The 
quantities A and uq, as well as the "working region" in the Borel parameter t are determined 
from the two-point sum rule for F 2 , and the sum rules for K and S do not contain free 
parameters. In practice, it turns out to be convenient to consider the ratio of three-point and 
two-point sum rules, which suppresses significantly all the sources of uncertainties. Adding 
a factor 4 (cf. Q2.24Q ) and taking into account the scale dependence, we obtain the following 
sum rule for the mass splitting between vector and pseudoscalar mesons: 



1 In this paper we conform to the standard duality region, which is a square in the (s, s') plane. Its 
size is fixed by the continuum threshold in the corresponding two-point sum rule. Our sum rules 
are not sensitive to the detailed shape of the duality region near the diagonal, which is different 
from the situation considered in |2(| . 



S 



a s (2t) J 



ds \df/e-WWpv(8, s') 
dse~ s/t p u {s) 



ml _ ml = S ( O-tSg f * . (3.13) 



o 



In the spectral functions, the renormalization scale is \x = It. A similar ratio determines the 
kinetic energy K of the heavy quark: 



Wq Wo 



2 dsJds'e~( s+s 'yWp K {s,s') 
K= ° ° 55 • (3.14) 



Jdse-^pu(s) 



Note that the explicit dependence on A has cancelled out. It is, however, present implicitly 
since the values of u and A are correlated. 

In the numerical calculations we use the following standard values of the vacuum con- 
densates (at the normalization point 1 GeV): 

(qq) = -(240 MeV) 3 , 



G ) = 0.012 GeV , 

7T / 

(qagGq) = (0.8 GeV 2 ) (qq) . (3.15) 

In this section we use the one-loop expression for the running coupling with Aqcd = 150 MeV 
and four active flavors. The corresponding values of the coupling are a s (lGeV) = 0.40 and 
a s (m B ) = 0.21. 

The results are shown in Fig. ||. For definiteness, we give values of the decay constant 
and the chromomagnetic splitting normalized at the scale of the B-meson mass. The sum 
rule for the two-point function is most stable for the value A = (0.4 — 0.5) GeV and ujq = 
(1.0 — 1.2) GeV, see Fig. |5|(a). The stability region starts already at values of the Borel 
parameter of order 0.3 GeV, and stretches practically to t — > oo. It is known, however, that 
stability at large values of the Borel parameter is not informative, since in this region the 
sum rule is very strongly affected by the continuum model. The usual criterium that both 
the higher order power corrections and the contribution of the continuum should not be very 
large (say, less than (30-50)%), restricts the working region considerably. In the particular 
case of the decay constant F(pj) one usually chooses 0.3 GeV < t < 1 GeV 0. In the case of 
the three-point functions, the sum rules are especially strongly affected by the subtraction 
of the continuum owing to the high dimension of the spectral densities. Thus, it is especially 
difficult in this case to ensure that the continuum contribution is sufficiently small, and this 
requirement forces one to work in a rather narrow region of the Borel parameter, close to the 
lowest possible value t ~ 0.3 GeV. This strategy is backed up by considerable experience of 



QCD sum rule calculations of higher-twist operators built of light quarks, see e.g. pl|j22]] . In 



this paper we take the working window in the Borel parameter for the three-point functions 
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( [3.1D and ( p.2|) to be 0.3 GeV < t < 0.6 GeV. We emphasize that this region should be fixed 
by considering the two-point function, and does not necessarily coincide with the stability 
plateau for three-point sum rules. The remaining instability should be considered as a part 
of the errors involved in the calculation. 

In Fig. |5](b) we plot the right-hand side of the sum rule for the heavy quark kinetic 
energy, Eq. (|3.14j ), and in Fig. |(c) the right-hand side of the sum rule (|3.12|) , as a function 



of the Borel parameter for two different values of the continuum threshold luq = 1.0 GeV 
and co> = 1.2 GeV. The sensitivity of the sum rules to the change of u) Q and of the Borel 
parameter (within the working region which is shown as shaded area) gives a conservative 
estimate for the accuracy of the results. We end up with the values 

{m 2 v - m 2 P )^= m ^ = (0.46 ± 0.14) GeV 2 , (3.16) 

K^ Q) = -(0.54 ±0.09) GeV 2 . (3.17) 

The given value of the mass splitting between vector and pseudoscalar mesons is our final 
result whereas the analysis of K will be extended to next-to-leading order in the next section 
(the superscript "LO" stands for "leading-order" result). The experimental values for the 
mass splitting in beautiful and charmed mesons are |23| 

m%. - m 2 B ~ 0.48 GeV 2 , 

m 2 D , - m 2 D ~ 0.56 GeV 2 . (3.18) 

These values are close to each other which indicates the smallness of higher order corrections. 
The measured splittings agree quite well with our result in (|3.16|) . 



IV. THE ORDER a s CORRECTIONS TO THE KINETIC ENERGY 

In this section we calculate the 0(a s ) corrections to the sum rule for the kinetic energy. 
This calculation is laborious, but worthwhile, since in the case of the coupling F the radiative 
corrections proved to be very large, see PJ24|. Following ||, we introduce the heavy-light 



current Jr, which is renormalization group invariant to two-loop accuracy (cf. (|2.14|) ), 

J r = (A)Wa s W-^« (l - S ?M\ , (4.1) 



with 

* = ^(^-^U-0.23. (4.2) 
The corresponding invariant coupling F is defined as (cf. fl2.12|) ): 



F(/i)a s (/i)-^W U-S^-Y 



(4.3) 
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The leading- and next-to-leading order anomalous dimensions of Jr, 70 and 71, as well as 
the coefficients of the /3-function are given in (|2.15 ). 



The two-point correlation function of the invariant currents Jr can be expressed in terms 
of invariant quantities and the strong coupling a s (fi) at the scale of the external momentum, 
\i = —2uj. A calculation of radiative corrections to the two-point function ( |3.4|) leads to the 
following sum rule for the invariant coupling F 0: 

F 2 e~ A/t = 



o 7 

4«s(2t)- W/3 ° / dss 2 e 



I 7T V 3 9 2 



- {l + ^ {2 - } 3 + ^ a s (2t)^-^^) ; (4 . 4) 

which is valid to two-loop accuracy. Here we have introduced the scale-invariant condensates 

.,(3) fJ?) 



a. 



7T 



C 4 = ( — G ) [l + 0(a.)] 



5 = (qtrgGqMa.M-^'W [1 + 0(a.)] . (4.5) 

The leading-order anomalous dimensions are 7/, = 270 = —8, 7q 5 ' ) = —4/3, the next-to- 
leading order anomalous dimension of the quark condensate equals = —404/9 + 40n//9. 
As a shorthand, we use A71 = 27! —71 = 704/9 — 1127r 2 /27. Apart from an overall scaling 
factor 2a s (2t)- 7 °/( 2/3 °), the sum rule for F in ( |4.4| ) differs from the one given in the last 
section by terms of 0(a s ). As for the determination of £ and K to one-loop accuracy, we 
neglected these terms and consistently used the leading-order spectral densities (|3.10|) in 
the sum rules ( |3.13| ) and ( |3.14| ). 



Our aim now is to calculate the double spectral function p K defined in Eq. (|3.8| ) to next- 
to-leading order. The calculation of the perturbative two-loop corrections to the diagram 
Fig. §(a) is the most difficult task. The necessary techniques are explained in App. [B], while 
the contributions of the individual diagrams in Fig. ^| are collected in App. |C|. The result 
reads: 

pert/ n / 3 \ 4 , f a s (/i) /41 4 2 7 2s\\ 



An 2 J ' { tt \ 9 9 2 fi 

a s {p) 1 



( s + s ')^2 + s /2_| s 2_ s /2|j_ (46) 



7T 27T 2 

In Feynman gauge, the spectral densities of individual diagrams contain terms like 
d/ds \n\s — s'\ which remind of the cancellation of infra-red divergences between partic- 
ular discontinuities, corresponding to real and virtual gluon emission. In the sum given 
above, however, they cancel. 
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The contribution of the quark condensate is easier to calculate. The relevant diagrams 
are collected in Fig. || and their spectral density is 



1 Os 
3 7T 



(qq) {s5(s') + s'S(s) + (s + s')5(s - s')} . 



(4.7) 



Finally, we take into account the contribution of the gluon condensate, see the diagrams in 
Fig. [|. The corresponding contribution to the spectral density equals 



(4.8) 



The leading tree-level contribution of the mixed condensate is given in (|3.10|) . We neglect the 
contribution of the four-quark condensate which is tiny in the case of F and of semileptonic 



form factors 25 



The results given above correspond to the calculation of the Feynman diagrams shown 
in Figs. 0, 3, and 4. To obtain the correlation function involving invariant currents, one has 
to multiply the spectral densities Eqs. (|4~6"1), (|4.7|) , and ( |4.8| ) by the corresponding power 
of the coupling. Then, after the Borel transformation, renormalization group improvement 
reduces to the substitution of the scale \x by 2t, twice the Borel parameter (see || for details). 
Combining all terms, we end up with the sum rule 



F 2 Ke~ K l t 



^ 



° S *S(8-S r ) 1 + 



a s (2t) /41 4 



7T- 



7T 



9 9 



-25 + ^ln* 
2 t 



a s (2t) 2 



7T 7T 



4 (s + s') (s 2 + s' 2 



I 2 12 \ 
\s — S 



3 71 



x I s5(s') + s'5(s) + (s + s')8(s - s ') 1 - 1 a s (2t)-^' Po O a 5(s - s') 



+ -a s (2t) 



( 7 (5) -2 7 o)/(2/3 ) 



5 5(s)5(s') 



(4.9) 



In both the sum rules ( |4.4|) and (|4.9|) the continuum contribution is subtracted also from the 
imaginary part of the running coupling a s (— 2uj). The term lns/t comes from the expansion 
of {a s {2s)/a s {2t)}-^/P° to first order §. 

In the numerical evaluation of ( |4.9| ) we use the values of the condensates as given in the 
last section and the two-loop formula for the running coupling, 



oi s (jj) 



An 



2/3 ln(/VAgf ) ; 



] _ ft ln(21n(/x/A^ ) )) 



Po 21n( / ,/AS| ) 



(4.10) 



with A^l = 260 MeV for four running quark flavors f23fl . As in the analysis of the leading- 



MS 



order sum rule ( |3.14[ ), we vary the Borel parameter t in the range 0.3 GeV < t < 0.6 GeV 
and the continuum threshold in the range 1 GeV < uj < 1.2 GeV. In Fig. |6|(a) we show 
the coupling F 2 , calculated according to flO]) and scaled up to fi = m B according to 
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(|4.3|). Figure ||(b) contains the kinetic energy K, calculated by taking the ratio of ( |4.9| ) and 
( |4.4|) . Both quantities are plotted as functions of the Borel parameter t for ujq = 1 GeV 
(A = 0.4 GeV) (solid lines) and u = 1.2 GeV (A = 0.5 GeV) (dashed lines). The working 
region is indicated by the shaded areas. As can be seen when comparing the two figures 
with each other, K is less sensitive to the value of ujq than F 2 , and we get the final result 

K = -(0.60 ± 0.10) GeV 2 . (4.11) 



The errors include the dependence of the sum rule C |4.9f) on variations of the Borel pa- 
rameter t and the continuum threshold uio. in the range of parameters considered, the 
non-perturbative terms contribute approximately 25%. 

The value ( |4.11| ) is very close to the leading order result (|3.17 ). The two-loop radiative 



corrections in the sum rules are large and change the value of the coupling F (ttib) by a 
factor two, cf. Figs. |5](a) and §|(a). Nevertheless, their net effect on the value of K is small 
due to a strong cancellation in the ratio of ( |4.9|) and fl4.4j ). It is interesting to note that the 
Coulombic radiative corrections, which contain an extra factor 7r 2 , are identical for the two- 
and the three-point sum rules. A similar effect was observed in the study of the Isgur-Wise 
function in PBI . 



V. DISCUSSION 

We have derived sum rules for three-point correlation functions in the heavy quark effec- 
tive theory, from which we obtain estimates for the matrix elements of the chromomagnetic 
interaction and the kinetic energy operator over heavy-light mesons. Our final results read 

(M(v)\K\M(v)) , N ^ o 

11 ' -(0.60 ±0.10) GeV 2 



(M(v)\h v h v \M(v)) 

(M(v)\S\M(v)Y =mB _ d 
(M(v)\h v h v \M(v)) ~ T 



-[(0.46 ±0.14) GeV 2 ] . (5.1) 



Note that the value of kinetic energy is renormalization group invariant, while the chro- 
momagnetic interaction depends on the heavy quark mass scale. The latter dependence is 
actually significant. The given value (in square brackets) at \x = rriB — 5.3 GeV corresponds 
to (0.60 ± 0.18) GeV 2 at the lower scale // = 1 GeV. 

Our result for the kinetic energy is two times lower than the value obtained in M from 
the expansion of two-point sum rules in QCD in the heavy quark limit. The approach of the 
present paper should be much more accurate, since the consideration of three-point functions 
allows one to suppress contaminating contributions of non-diagonal transitions. Still, this 
disagreement is disturbing, since QCD sum rules normally have a (20-30)% accuracy. The 
large value of the kinetic energy obtained in Q may be an artifact of the used procedure 
with a redefinition of the Borel parameter by terms of order l/m,Q. From our point of view, 
no such redefinitions are allowed. 

The corrections to B-meson masses, corresponding to (|5.1|) , are 

m B = m b + A + —[(0.25 ± 0.20) GeV 2 ] , 
2rrih 
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m B , = m b + A + [(0.70 ± 0.20) GeV 2 

2m h 



(5.2) 



respectively. Note that the 1/rrib correction to the pseudoscalar B-meson is very small. 
This gives strong support to the procedure of the evaluation of the 1/rriQ correction to the 
leptonic decay constant in ||, where the 1/wiq correction to the pseudoscalar meson mass 
have been neglected. For completeness, we quote the result of M 



F(m Q ) 



1 



(0.8 - l.l)GeV 

m Q 



(5.3) 



which agrees to all other QCD sum rule calulations f2"7| , |2"8"|] , with the only exception of ||, 
where a much larger correction is claimed. 

The phenomenological applications of our results include l/rnq corrections to semilep- 
tonic form factors and inclusive B-decays. The rather small value of the kinetic energy 
obtained in this paper may require the reconsideration of recent estimates of these correc- 
tions in pi . 
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APPENDIX A: THE l/m Q EXPANSION OF CORRELATION FUNCTIONS 

The aim of this Appendix is to establish a formal connection between the I/tuq expansion 
of correlation functions in QCD, and the corresponding correlation functions in the HQET. 
Let us consider the following two-point correlation functions of vector and pseudoscalar 
currents: 



IW<?) = i \ rf 4 xe^(0|T{g(x)7 M ^)K0)7^(0)} |0> 



-9^ + ^ K(g 2 ) + M^n s (g 2 ) 



n 5 (g 2 ) = % J d 4 xe iC!X (0\T{q(x)z l5 b(x)b(0)i l5 q(0)} |0> . (Al) 

The correlation functions n„(g 2 ) and n 5 (g 2 ) have poles at q 2 = my and q 2 = nip, respec- 
tively, and the contribution of the ground state mesons in the vector chanel, V, and in the 
pseudoscalar chanel, P, equals 



Ii v {q 2 ) = 4^- 2+ ... 
m v — q l 
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where the couplings fp and fy are defined as in Q2.16Q . Our objective is to make a systematic 
expansion of the correlation functions in ( |A"l] ) near the poles (|A^) in powers of the large 
quark mass mq. Each correlation function has symmetric pairs of poles at go — i^M, 
which correspond in an obvious way to contributions of particles and antiparticles. For 
definiteness, we consider the expansion near the particle-type discontinuity, and to this end 
define a new variable uj by 

q,t = v li (mQ+u)). (A3) 
The physical couplings and masses of the mesons are expanded as 



fM\pUM — Fm + SFm, 

m M = m Q + A + 5m M , (A4) 

where 5Fm and Stum are of order I/tuq. We anticipate here that the mass splitting between 
pseudoscalar, M = P, and vector mesons, M — V, is an I/tuq effect. Note that as a quite 
general feature of QCD, the extraction of asymptotic behavior (in our case in the heavy 
quark mass) introduces divergences. The separation of the asymptotic value of the coupling 
F tacitly assumes some regularization and the summation of logarithms of the heavy quark 
mass, which is made explicit using the formalism of the heavy quark effective theory. 

Putting together Eqs. (|A2|) , (|A3[) , and (|A4j) , we obtain, e.g. for the vector correlation 
function: 

/LmL F 2 25mF 2 2F5F , 

H h non-singular terms, (A5) 



m 2 B , - q 2 2A (2A) 2 2A 



where 



A = A - u. (A6) 



The expansion of the correlation functions in ( |Al| ) is not so immediate. Using the heavy 
quark expansion of currents, Eq. ( [2.11| ), and the 1/mq expansion of the Lagrangian, 
Eq. (|2.3|) , we obtain for the vector correlation function 

d 4 xe^ m ^\0\T{q(x) 1 ,Q(x)Q(0) lu q(0)} |0> = 

i J rf 4 xe-^(0|T{g(x) 7 >,(x)^(0) 7 ^(0)} \0)0^) 

l —i J rf 4 xe^(0|T{g(x) 7 ^(^ ± )^(x)^(0) 7 ^(0)}|0)^= m ^ 



3 



+ J rf 4 xe i ^(0|T|g(x)7>,(x)^(0)(-^ ± ) 7 ^(0)||0)^= m ^ 

+ d 4 x d 4 ye <w ^(0|T{g(x) 7 ^(x)/C(y)fc t ,(0) 7 j-g(0)} \0)^= m ^ 



15 



2m c 



—i 2 f rf 4 xAe M ^ ) (0|T{g(x)7>,(x) l S( Z /)^(0) 7 ^(0)}|0)^= m ^ 



« s (wq) 



+ longitudinal part + O + o(l/wf ). (A7) 

* 2mn 7T ' v 



Most of the correlation functions in ( |A7|) have been introduced already in the main text. In 
addition to them, we encounter one more invariant function, 

i J d 4 a;e <WW! (0|r{g(a;)r 1 (iDi)^( a ;)^(0)r 2g (0)} |0) = -^{r^r^n^u;), (A8) 

which can be related, however, to the invariant function H(u) in (|3.4j ), choosing 1^ = j a , 
integrating by parts, and using the equations of motion for the effective heavy quark field 

(vD)h + = and for the light antiquark q{—i ft) = 0. One obtains Q: 

H 1 (lj) = ~uII(u). (A9) 
Combining everything we arrive at 

n,(g 2 = (mg + uof) = (l- ^^9l) U (cu) + -^U(u) + —T k (uj, co) - — T s (uj, u). 

V 3 7r J 3m Q m Q m Q 

(A10) 

The expansion of the correlation function of pseudoscalar currents is obtained by similar 
manipulations. The result reads 

n 5 (g 2 = (mg + uf) = ( i - tS!^9l) n( w ) + — n(«) + — r^(^, w ) + — t s (^, w ). 

(All) 

Note that these relations are valid up to constant terms only, which come from distances 
~ 1/ mg and are lost in the heavy quark expansion. Using the operator product expansion 
results for the vector and pseudoscalar correlation functions in QCD, it is easy to verify that 
the answers for the contributions to the correlation functions in HQET of the quark and the 
mixed quark-gluon condensate given in the main text indeed satisfy the general relations in 
( |A11| ) and ( |A10| ). For contributions of the pertubation theory and the gluon condensate this 
check is not so immediate, since in the process of calculation we have discarded contributions 
which vanish after the double Borel transformation. 

The last step is to insert a complete set of meson states in the correlation functions 
and to separate the contribution of the lowest energy. Note that at uo = uo' the three-point 
correlation functions Tk(u),u') and Tz(u,u) r ) contain contributions with both a double pole 
and a single pole at uj = A: 

F 2 {n)K 2F 2 (fi)G K . , 

T K (u, u) = j— j \Y2 + -^71 r + non-singular terms, 

|2(A — u)\ z 2{A — uj) 

T s (cj, u) = t—j \Y2 + ~o71 T + non-singular terms, (A12) 

[2{A — uj)\ z 2{A — uj) 
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where K and E are matrix elements of the kinetic energy and the chromomagnetic interaction 
operator over the heavy-light meson in the effective theory. The quantities Gk and are 
defined as 

(0\iJdx]C{x)h v {0)Tq{0)\M{v)) = ^F{/i)G K {ii)Tr{TM{v)} , 

(0| % J dxS(x)h v (0)Tq(0)\M(v)) = l -F(ii)G^)2d M Ti{YM(v)} , (A13) 

where it is implied that the matrix elmements of the two-point correlation functions on the 
left-hand side include summation over all intermediate states except for the lowest-lying 
one. For example, 

(0|i / dx IC(x)h v (0)Tq(0)\M(v)} = 

= J2 m\M'(v)) } - - {M'(v)\h v Tq\M(v)) (A14) 
m>+m Zm B (A -A) 

where A' = itlm' — ttlq. Up to this remark, our formulae coincide with the ones derived in 
I- 

Collecting everything and comparing to the expansion of the relativistic expression in 
Eq. flA|), we obtain 

5m = -— [K + ef M £] , 

SF _ A 1 d M 

F 6m Q 2m Q m Q 

2mp(mv — mp) = my — mp = 4E, (A15) 
which is the desired result. The relations in (|A15|) can of course be derived in a variety of 



ways and are by no means new. Our derivation is suited best for the application of the sum 
rule technique to the evaluation of relevant parameters. 

APPENDIX B: CALCULATION OF SPECTRAL DENSITIES TO TWO LOOP 

ACCURACY 

In this appendix we give all the formulae necessary to calculate the spectral densities 
of the correlation functions (|Q| ) and ( |3.2|) to two-loop accuracy. To this end, we use 
a technique proposed by |30| which amounts to calculate double spectral densities by a 



repeated application of the Borel transformation to the amplitude itself. Let A be some 
amplitude depending on the virtualities uj and uj' , which we want to write as a dispersion 
relation (subtracted if necessary) in both variables: 

A(uj, uj') = fds ds' P \. , + subtractions. (Bl) 

J (s — uj)(s' — UJ) 
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Let us assume that the amplitude A(u,u') is calculated, e.g. in the from of a Feynman 
parameter integral. Then, applying a Borel transformation in uj and cu', we obtain, as the 
second step: 

4 (*) 4/(0^(^0 = ^(M') = fdsds' p(s,s')e- s/t - s,/t '. (B2) 
Finally, we apply a Borel transformation in r = 1/t and r' — l/t' to get 

4 (1/s) 4' i (1/r, 1/r') = ss'p(s, s'), (B3) 

which is (apart from a trivial factor) just the desired spectral density. 
Here we made repeated use of the formula (with w B = — uj) 

4(i)e~ Qa,B = 5(1 -at). (B4) 

We spend one subsection to derive necessary formulae for each of the above three steps, and 
finally illustrate the whole procedure by the sample calculation of a one and a two-loop 
diagram. 



1. Reduction of Loop Integrals to One Parameter Integrals 

The general one-loop integral is defined as 

r d D k 1 
h(a,p,q) = y^-p [k 2 ] a (u + vky(u' + vky ^ 

Using the parametrisation 

1 - T{a + P) %^^ (B6) 



A a PP T(a)T(p)J (A + yP) a +P 

we find 



h(a, p, q) = _1_ j( a> p> q) F{q _ i ; p _ i ; 2a + p + q - D) (B7) 



with 



Tf \ T{2a + p + q-D)V{D/2-a) 

J «, P, Q) = r( ypv ypv x , B8 

r(a)r(p)r(g) 

i 

F(a, 6, c) = Jdz z a (1 - zf [(u' E z + uj e (1 - z)}~ c , (B9) 
o 

where we have already continued to Euclidean space: uje = —uj, oj' e = — uj' . In the limit 
q — > we find 
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with 



J 1 (a,p l O) = (-ir +P i ^/(a,p)^- 2 ^ P (BIO) 



T( \ T(2a+p-D)T(D/2-a) 

J(a ' p) = fWlp) (B11) 



and a corresponding formula for p — ► 0. 

The general two-loop integral can be expressed as 

I 2 (a,fl,7,p,q,r,s) = 

dPk d D l 1 



(B12) 



(2tt) d (2tt) d [k 2 \ a [l 2 f [{k - I) 2 }"' {u + vk)P {uj + vl)i (a/ + vk) r (a/ + viy 
Simple cases are 

h(oc, P,0,p, q, r,s) = h(a,p, r) h(/3, q, s), (B13) 

I 1 \a+"/+p+r+s+l 

7 2 (a,0,7,p,0,r,s) = 4a+ln D J (7,s) J{j,p,r + s + 2a - D) 

xF(r + s + 2a - D - l,p - l,2(a + j- D) +p + r + s), (B14) 

r, n n n n gg/2 - a)T(D/2 - 7 )r(2(q + 7 - g) + g + g) 
J 2 (a, 0, 7, 0, g, 0, s) = — 

r(a)r(7)r(g)r(s) 



40+771-° 



F(s-l,g-l,2(a + 7-D) + g + s), (B15) 



I 2 {a, 0, 7, P, 0, 0, s) = ( l) *^D +l J (7> s) J(a, P, s + 2 7 - £>) 

xF(s + 27 - L> - l,p - 1, 2(a + 7 - D) +p + s), (B16) 

f_l N ) Q +/ 3 +7+P+''+l 

J 2 (a,/3,7,P,0,r, 0) = 1 Aa+0+ ^ D F(r - l,p - 1, 2(a + /3 + 7 - L>) +p + r) 

s r(fl + 7 - D/2)T(D/2 - (3)T(D/2 - 7) 

Similar formulae can be obtained be exchanging u; and a/ or the loop momenta. The special 
case r = s = is considered in |~7| where recursion relations for I 2 were derived by using 
the method of integration by parts (cf. |]l7lj3l|). We employ the same technique to relate 
the general I 2 to the special cases solved above and find the following recursion relations: 

I 2 (a,0,7,p,q,r,s) = 

= o — , 1 n { 1 /2 ( a ' °> 7>P - 1> 9 + 1> ^ «) + « J 2(«, 0, 7, p, g, r - 1, s + 1) }, (B18) 

27 + g + s - D 
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J 2 («,0,7,0,g,r, s) = - — { qh(a, 0, 7, 0, q + 1, r - l,s) 

^7 + s — i-/ 

- g/ 2 (a, 0, 7, 0, g + 1, r, s - 1) + s I 2 (a, 0, 7, 0, g, r - 1, s + 1) }, (B19) 
I 2 (a, (3, 7, 0, g, r, 0) = ^— ^ ■ — {(3 1 2 (a, (3 + 1, 7 - 1, 0, g, r, 0) 

^ - 7 - u (2/9 + 7 + g-£>) 

-(3I 2 (a- 1,(3 + 1,7,0, g,r, 0) - g/ 2 (a, /3, 7, 0, q + 1, r - 1, 0) 
+ 7(1-^) [/ 2 («-l,A7 + l,0,g,r,0)-/ 2 (a,/5-l, 7 + l,0,g,r,0)]}, 



(B20) 



I 2 (a,P,-f,p,q,r,0) = — — 1 {/3 J 2 (a, /3 + 1, 7 - l,p,g,r,0) 

L> — Z7 — p — g 

- /? /.(a - 1, /3 + 1, 7 ,p, g, r, 0) - g J 2 (a, /3, 7 ,p - 1, g + 1, r, 0)}. (B21) 

Note that one has I 2 (a, (3, j,p, g, r, s) = I 2 (/3, a, 7, g,p, s, r). 
Finally, we encounter the integrals 

f M e 

7 (2tt) 2D l 2 {l-k) 2 {uj + vk){uj' + vk){uj' + viy 1 ; 

/ ^fcrf^ k 2 

7 (2tt) 2 ^ Z 2 (Z- kf{uj + vk){uj + vl){uo' + vk){u' + viy K ' 

After a change in variables, k — > k + 1, and noting that the integral over I can depend on v 
only, we find 

J {2ir) 2D l 2 k 2 {uj> + vk){uj + vk + vl){u' + vk + vl) ^ ' 

xuj E h(2- D/2, 1,1). (B25) 
Similarly, one obtains 

£? = 2cu E {u' E J 2 (0, 1, 1, 1, 1, 0, 1) - J 2 (0, 1, 1, 1, 0, 1, 0) - cu E / 2 (0, 1, 1, 1, 1, 1, 0)}. (B26) 



2. First Application of the Borel Transformation 

In this appendix we are concerned with the calculation of 

K(p, q; a, (3, 7) = B u (t) 4/(0 uo E u'i F(p, q, 7). (B27) 
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For non-negative integer a, (3 we can write 

Kip, q, a, (3, 7) = B u (t) B^(t')^u% u% jdz z p (1 - zf Jds s^ 1 e -(^^(i-)) 



f_V\a+l3 jq+/3 

t a t' p — t 1+a+p -^ t n+ P +q ^(t + t'y-P-i- 2 . (B28) 



r( 7 ) dt a dt*p 

Likewise we need K(0, 0, a N, (3 e N, 7). In the case 1 + /3 + [a] < 0, where the Gaussian 
bracket [a] denotes the integer part of a, one finds 

K (0, 0, a £ N, /3 e N, 7) 

* ^ 4 77777^7 / <** ^ (* " r^) • (B29) 



r(-«)r(7) rft'/ 3 + •/ v t + t' 

t+t' 

The integral apparently becomes singular for a > 0. It is convenient to collect divergent 
terms in T functions rather than in terms of hypergeometric functions. This can be achieved 
performing (1 + (3 + [a]) integrations by parts, if this number is greater than zero. Defining 



1 

A(a,(3) = J dzz a ( 



z 



t + t' 



t+t' 

1 t a t'? +1 



2 F 1 (-a,l + (3,2 + (3,-t'/t) (B30) 



l + (3 (t + t') 1+a +P 
we find, if P denotes an operator for partial integration: 

P A(a, p) = T ^(t + t')- {l+P) t'^ 1 - ^ A(a - 1, + 1), 

le A{a ' p) = p Jt+Hy A{a > (B31) 

This yields finally (for (1 + (3 + [a]) > 0) 

{-If 



K(0,0,a£N,(3 e N, 7 ) 



r(-a)r( 7 ) 

,0 ,/l +a+i a_ 7 J, . -(l+a) 

x ^ 77« 7 777" [P] 1+/?+H / rf^ 1+Q - 7 U - ) • (B32) 

dt'P (t + t') 1+a L J 7 V t + t' J 



t 

t+t' 



21 



3. Second Application of the Borel Transformation 

In terms of the new variables r = l/t and r' — 1/t', a typical contribution to the Borel 
transformed amplitude is of the generic form 

1 (B33) 



with integer a, /3, 7. Some terms are accompanied by logarithms of r, r' or r + r', which 
we take into account by replacing a by a + e, e.g., and expanding around e = afterwards. 
The expressions we are dealing with have mass dimension five, thus it is feasible by repeated 
replacements r — > (r + r') — t' to reduce all powers in the denominator to non-negative 
values (plus e for logarithms). We are left with the calculation of (with in general non-integer 
values of a, j3, 7) 

L(a, (3, 7) = B T (l/s) BAl/s') + r , )7 • (B34) 

As usual, we express the denominator in terms of integrals of exponential functions that 
have a simple behaviour under the Borel transformation (cf. Eq. (|B4|) ) : 

L(a,/3,7) = 

00 

= B T (1/S) B T/ (1/S') I d Ul du 2 du 3 e-u^r'-u^r+r') 

r(a)r(/3)r(7) j 

e(s' -s) J duu^ 1 (s - u) a ~ l (s' - uf- 1 + (s^s',a^ (3) 



s s' 



i»r(/3)r( 7 ) 







1 f s «+7 8 'f> J» 2j p x (i _ % a + 7> s / s ') + ( s <_> s / j a <_> . (B35) 



T(a)T(0) { I> + 7 ) 

As special cases we find 



ss' 



B T {l/s) BAl/s') - = - °1 0(« - (5 - s'T- 1 s'^\ (B36) 

r a (T + r) 7 r(ajr(7j 

~ 1 <?<?' 

S T (l/ a ) S T ,(1/ S ') -— ^ = ^ S ^ *(* - *')■ (B37) 

Fortunately enough, in the sum of all diagrams most logarithmic terms cancel and the 
remaining ones are of type 

ln(r + r'), (B38) 
so we need not expand the hypergeometric function in (|B35|) . 
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4. A Sample Calculation 

In this appendix we calculate the diagrams D a , Fig. ||(a), and D b , Fig. 0(b), which 
contribute to the kinetic energy. We use dimensional regularization in D = 4 + 2e dimensions 



and work in Feynman gauge. We use T\ = T2 = ^75 (cf. Eq. (|3.1|) ) and do the traces 
explicitly, which adds a factor two. Thus, each diagram contributes to the spectral density 
px, Eq. (|3.14j ), with a weight factor 1/2. 



In momentum-space, the lowest order triangle diagram D a can be expressed as 

a ~ J(2tt) d k*(u + vk)(u' + vk)' 1 j 

Since we are interested in the Borel transformed expression only, we write vk = (vk + u) —uo 
and keep terms with non-polynomial dependence on both u and u'\ 

D a = -MBUt)B w ,(t') f4%{" 2 , T J, FT ZS +^ TTf T L FT h X ( B4 °) 

J (2,71) y (lj + vk){u' + vk) k z (uj + vk)\u' + vk) J 

The first loop-integral vanishes, thus we are left with 

b. - - «*.(«> *,(*V /<g> ^7T^7T^) - (B41) 

The integral encountered here is a special case of the general one-loop integral 

1 

^(1,1,1) = —^^ Jdx(u' E x + u E (l-x)) D -\ (B42) 



where we have changed to the Euclidean variables ue = —oj, io' E = —oj'. As a shorthand, 
we use the notations g\ = T((D — 4)/2) and gi = (4 — D)T(A — D). Although the integral 
can easily be solved in the present case, the above form is convenient for applying the Borel 
transformation afterwards. In general, we have to apply the Borel transformation to terms 
like 

1 

K(j>,q;a,P,i) = B u {t) B w ,{t') u) a E u)% Jdxx p (1 - x) q {u' E x + u E {l - x))-\ (B43) 



With the formulae given in App. [TJ2 we find 

K(0, 0; 3, 0, 4 — D) = — (D - A)(D - 2)(D - l)Dt D t' D (t + t'Y 1 - . (B44) 

92 



Taking altogether, we have 

3 

47^/ 2 

5 

I +0(e). (B45) 



D a = - 91 (D - 4)p - 2){D - l)D t D t' D (t + if)-^ D 



36 1 / tt' 

T^tP [t + t> 
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Finally, to calculate the double spectral density of D a , we apply to D a a Borel transformation 
once again. More accurately: if p a is the spectral density, we have 

Pa(s, s') = — B T {l/s) B T ,(l/s') \t) a [s -> 1/r, s' -> 1/r') 

= -4^ 9 'r^) sD ^- s '»- (B46) 

We have checked explicitly that in D = 4 dimensions this result agrees with that obtained 
by Cutkosky rules. 

Next we turn to the two-loop diagram given by 

D b = B u (t) 4/(0 327ra s { 2u E u' E {u E + lo' e ) I 2 (l, 1,1,0, 1, 1, 0) 

- l - (3uj e + u' E ) I 2 (0, 1, 1, 0, 1, 1,0) - | (3a4 + J 2 (l, 0, 1, 0, 1, 1, 0) 

+ ^ + o4)/ 2 (l, 1,0,0, 1,1,0)} (B47) 
with the two-loop integrals 

' ' (B48) 



(2tt) d (2tt) d [fc 2 ] Q [Z 2 ]/ 3 [(Jfe - l) 2 p (lo + vk)P {uj + vl)i (a/ + vk) r (u/ + vl) s ' 

In the above case all integrals can be related to one-loop integrals, except for 
7 2 (1, 1, 1, 0, 1, 1, 0). Here it proves useful to employ the method of integration by parts 
From 

fJ^J^L± = (B49) 

J {2n) D (2tt) d dip [k 2 } a [l 2 Y[(k - l) 2 }^{u + vl)"{uj' + vk) r 1 ' 

we get 

( J D-2/?-7-g)/ 2 («,/3, 7 ,0,g,r, 0) - 7 J 2 (a, /3 - 1, 7 + 1, 0, q, r, 0) 

+ 7 I 2 {a - 1, 0, 7 + 1, 0, g, r) + gcu J 2 (a, /3, 7 , 0, g + 1, r, 0) = 0. (B50) 

This relation still contains the term i 2 (a, j3, 7 , 0, q + 1, r, 0), which is as difficult to calculate 
as the original expression. We can, however, make use of a corresponding relation that 
follows from 

[J^LJ^L± K = (B5i) 

J {2ti) d {2ti) d di^[k 2 } a [i 2 Y[{k-i) 2 ^{uj + viY{uj' + vky K ' 

to eliminate 7 2 (a, /3, 7 , 0, q + 1, r, 0) and get 



24 



I 2 (a,/3,j,0,q,r, 0) = - ^^r— — ■ — {(3 1 2 (a, (3 + l,j - 1,0, q,r, 0) 

P - 7 - ~ (2/9 + 7 + <? - £>) 

-Ph(a- + 7,0,g,r, 0) - g/ 2 (a, A 7 , 0, g + l,r - 1, 0) 

+ 7 fl - ^) [ h{<* - 1, A 7 + 1, 0, g, r, 0) - I 2 {a, (3 - 1, 7 + 1, 0, g, r, 0) ] }. (B52) 

Using this expression, one can relate Di, to K, and apply the second Borel transformation 
according to the formulae given in App. [B]3. Note that both the diagram and its double 
spectral function are finite in the limit D — > 4, whereas all the other two-loop diagrams in 
Fig. |2| have to be renormalized. The calculation of the remaining diagrams proceeds along 
the same lines. The results are given in App. y. 



APPENDIX C: EXPRESSIONS FOR THE TWO LOOP DIAGRAMS 



In this appendix we collect the explicit expressions for the renormalized Feynman dia- 
grams in Fig. ^ after application of the Borel transformation in both external momenta. 
The missing diagrams with a gluon line connecting the operator vertex with a heavy quark 
line vanish in the limit of equal heavy quark velocities. D a , Dj,, . . . , denote the Borel trans- 
formed diagrams in Fig. 0(a), (b), . . . , respectively. Each diagram contributes with a weight 
1/2 to the spectral density Eq. (|4.6|) . 

br = W7)v { u ' (t + tr ~ 2(3 111 th (t ' 2 + m ' + m ' 2) 

-2(' 3 ln-^ 7 (i' 2 + 5i«' + 10« 2 )), (C2) 

= ^ J^y^ {-tt'(t + t')* - 12ft»(t + t'f 

+ 24t 3 t' 3 ( - - — + In 2 + In t -^—) + 2t 3 In (t 3 + Qt 2 t' + 15tt /2 + 22t' 3 ) 
\4 3 [i J t + t' y ' 

+ 2t' 3 In (t r3 + Qt'H + Iht't 2 + 22t 3 ) | , (C3) 

Dr = ^^^{-(t + to 4 + tnt + to 2 



+ 12tY 2 + lE - In 2 - In t±£\ } , (C4) 
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D p g cn = D p f CIt (t^t'), (C6) 
^ = ^^7^, f-l + 4 7E -41n2-41nV^V (C7) 

* 7T 7T 2 (t + t') 5 V pL(t + t') J V 7 
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FIG. 5. The sum rules ( p. 12 ) to leading-order accuracy (i.e. using the spectral densities ( |3.10[ )) 
as functions of the Borel parameter t for different values of the continuum threshold ujq (solid 
lines: u = 1.0 GeV [A = 0.4 GeV], dashed lines: w = 1.2 GeV [A = 0.5 GeV]). The shaded areas 
indicate the working regions of the sum rules, (a) F 2 (/j, = mg), (b) the kinetic energy K, (c) 
my — mp = 4S at the scale [x = tub- 
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FIG. 6. The sum rules (|4.4Q and ( |4.9|) in next-to-leading order accuracy. The parameter values 
are the same as in Fig. [|. (a) F 2 (pL = tub), (b) the kinetic energy K. 
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